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Abstract 

Ammunition deterioration during storage has considerable eco- 
nomic consequences. A reliable prediction model for the ammunition 
deterioration rate is necessary for long-term procurement and main- 
tenance planning. A random effect growth curve analysis is employed 
to formulate a prediction model for ammunition deterioration rates in 
terms of concurrent characteristics such as depot condition and vendor 
information. The resultant prediction model can be used to determine 
the appropriate time for reorder or renovation of ammunition before 
the quality reaches unacceptable levels. A two-stage analysis is used 
to estimate parameters involved in the prediction model. Necessary 
estimation methods are discussed. An example is given to illustrate 
the implementation procedure of the prediction model suggested in 
this paper. 

Key Words: Random Effect Logistic Regression Model, Deterioration 
Rate, A Two-Stage Estimation 



1 



1. INTRODUCTION 

Each year, the U.S. military services purchase millions of dollars worth 
of major caliber ammunition from several vendors to prepare for possible 
war. Ammunition lots purchased are stored in depots until they are used. In 
view of the fact that most ammunition is produced long before its ultimate 
consumption, it is important that the material be adequately stored and 
maintained to remain usable. 

However, according to the GAO report (The Comptroller General (1982)), 
the military services are experiencing significant problems in maintenance 
of depot and renovation of ammunition. These problems cause ammunition 
to deteriorate faster than necessary maintenance can be performed and a 
serious backlog of ammunition in need of maintenance. In an attempt to 
prevent such undesirable phenomena, in this paper, a quantitative model 
that can relate a depot maintenance plan to the quality control of ammu- 
nition is considered. 

The objective of this paper is to outline a prediction model for ammuni- 
tion deterioration in terms of depot characteristics and vendor information. 
The goal is to provide estimation methods by which one can determine the 
appropriate time for reorder or renovation of ammunition before the qual- 
ity reaches unacceptable levels. The suggested model can also be used for 
vendor control and examination of depot maintenance policy. 

In order to formulate such a prediction model for ammunition dete- 
rioration rates, a random effect growth curve analysis is employed. The 
analysis of a random effect growth curve model has been extensively used 
in biomedical research (Rao (1965), Stiratelli et al. (1984), and Vonesh 
and Carter (1987)). An application of growth curve analysis to physical 
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engineering is given in Sohn and Mazumdar (1991). In general, the goal of 
growth curve analysis is to estimate random growth rates of a time varying 
principal attribute in terms of concurrent individual characteristics which 
do not change over the experimental period. The analysis of growth curves 
is typically carried out by using a two-stage random effect model. The 
first stage consists of the within-individual regression model in which the 
within-individual growth rates for the serial observations are estimated for 
all individuals. In the second stage, estimated growth rates are related to 
a set of covariates representing concurrent individual characteristics. 

In this paper, a growing pattern of defective proportion of an attribute 
of ammunition lot over time (deterioration rate) is viewed in analogy to 
a growth rate of a principal attribute of an individual experimental unit 
while vendor information and depot characteristic;- associated with each 
lot can be considered as concurrent individual characteristics. In order to 
predict deterioration rates of a certain type of ammunit ion in terms of depot 
characteristics and vendor information, a random effect logistic regression 
model (Stiratelli et al. (1984)) is employed assuming that the pattern of 
deterioration follows a logistic function. For estimation of parameters in 
a random effect logistic regression model, a two-stage estimation based on 
the method used in Korn and Whittemore (1979) is employed. An example 
is given to illustrate the implementation procedure of the random effect 
logistic model suggested. 
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2. RELATED STUDIES 

The most important deterioration processes occurring to ammunition 
during storage are corrosion of the metallic components and chemical and 
physical changes within the internal structure of the explosives, which may 
result in reduced stability and performance (Eckman (1974), Eriksen and 
Strpmspe (1978), and Wilson and Bird (1986)). The factors which have the 
greatest influence on these processes are temperature and humidity which 
are directly related to the location and condition of the depot. Even under 
similar conditions, different types of ammunition may degrade with differ- 
ent deterioration rates. In some instances, different vendors provide the 
same types of ammunition that, nonetheless, have varying deterioration 
rates. Most of the studies in the literature do not reveal the quantita- 
tive relationship between deterioration rates and depot condition, types of 
ammunition and vendor information. 

For instance, Eriksen and Strpmspe (1978) conducted an extensive ex- 
periment covering a 9 year period to examine deterioration of six types of 
artillery ammunition (105 mm HE, 90 mm HE, 57 mm HE, 40 mm HE, 
hand grenades Mk2, and antitank mines M6A2) at three aboveground and 
three underground depot sites. Various attributes were observed: vacu- 
um stability and water content of propellant, impact sensitivity of primers, 
burning rate of propellant, muzzle velocity, and the number of misfires and 
duds. Inspection frequencies varied from two to four times covering com- 
binations of starting year (1968) and the subsequent 1st, 5th, 7th and 9th 
year after storage, depending upon the attributes inspected. Based on these 
repeated measurements, the contingency table and the analysis of variance 
were used to test if the levels of ammunition deterioration were significantly 



4 



different in terms of the conditions of depot and the types of ammunition. 

The results of their study are as follows: Overall, aboveground sites 
surrounded by inland climate performed better than underground depots 
that were old and technically in poor shape; the visual inspection revealed a 
considerable difference in corrosion attack between the ammunition types; 
water content both in the primer and the propellant charges varied signif- 
icantly depending on the depot in which lots of ammunition were stored; 
most of the chemical analyses were of little interest and should be postponed 
to the end of trials. Their approach, however, did not provide a compre- 
hensive model which related the level of deterioration of certain type of 
ammunition to depot condition. If the level of deterioration over time can 
be expressed as a function of depot conditions or/and vendor information, 
reorder time of a certain type of ammunition which reaches certain quality 
can be predicted in terms of given information at the depot. An appli- 
cation of growth curve analysis which facilitates such a prognostication is 
illustrated in the following section. 
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3. A GROWTH CURVE ANALYSIS FOR AMMUNITION 
DETERIORATION RATE 

Consider a certain type of ammunition (say, 105mm HE). In this section, 
the assumptions used in the growth curve analysis to predict the deteriora- 
tion rate of 105mm HE ammunition are summarized. Similar assumptions 
can be made for other types of ammunition. 

Basic unit of purchase is a lot and as a result of acceptance sampling, 
the qualities of incoming ammunition lots would be homogeneous regardless 
of different sources of manufacturers. While several hundreds of such lots 
are stored in various locations, the quality of ammunition deteriorates with 
the passage of time. Deterioration rates may vary depending upon both the 
condition of depot in which ammunition lots are stored and manufacturers 
who provide the ammunition lots. 

In order to inspect the pattern of deterioration of various lots over the 
time, it is assumed that a set of initially selected lots in various depot- 
s is marked and these lots are repeatedly inspected on a sampling basis 
without rectification over a certain time period. Lot size is assumed to 
be sufficiently large that it would not be depleted at the end of repeated 
sampling inspections. Modes of inspection can cover visual inspection, lab- 
oratory examination and functional trials (Eriksen and Strpmsoe (1978)). 
The proportion of samples which do not satisfy a predetermined level of an 
attribute or the combinations of attributes can be used as a proxy for the 
accumulated level of deterioration at the time of inspection. For instance, 
among many other attributes of functional trials, the relative number of 
misfires and duds to the sample size observed at each inspection cycle can 
be considered as a measure for the level of deterioration of each lot at the 
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time of inspection. 

For each lot, the proportion of defective rounds is assumed to increase 
as time goes on following a logistic curve. The logistic model is well es- 
tablished as a basis for analyzing such a phenomenon when the response 
variable follows binomial distribution. The logistic model formed to de- 
scribe the deteriorating pattern for each lot is called a within-individual 
model. Deterioration rate of a lot of ammunition, which represents a dete- 
riorating pattern over time of a lot often differs from that of another lot. 
The model in which varying deterioration rates between lots are compared 
in relation to depot characteristics and vendor information associated with 
each lot is called the between-individual regression model. Depot condition- 
s and vendor information associated with each lot of certain ammunition 
are assumed to remain constant during the experimental period while they 
would vary over different lots. Necessary notations for a two-stage random 
effect logistic model for ammunition deterioration rate is as follows: 
NOTATIONS 

v. index for individual lot, i = 1 , N 

N: total number of ammunition lots employed in the experiment 

j: index for inspection, j — 1, .., n t 

n t : the number of inspections to be done on a lot i over experimental 
period 

Uj : jth inspection time of a lot i 

Yij\ sample size of lot i at the jth inspection 

Xijk'. a dichotomous variable indicating if each item k taken from the 
sample at the yth inspection of lot i is defective or not 

Xif the number of defectives found in a sample size of X i j='£,k'Li X xjk 
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Pq = Xij/ Yij\ proportion defective at the jth inspection of a lot i 
exp(/?o)/(l+exp(/?o))‘- initial proportion defective of ammunition lot 
Pi : the deterioration rate of lot i (log odds ratio over the omit time 
interval for lot i) 

{.ZiiiZa, ■•jZim)'- covariates indicating concurrent characteristics of lot i 
(e.g., dummy variables indicating different combinations of depot condition 
and vendor information of ammunition lot i ) 

7b -1 7 m- regression coefficients for (z,i, z, m ) 

m: total number of covariates in the between-individual model 

T a pxl vector of 72 ’s, (71, 7 m )' 

Z: aniVxm matrix of covariates ( z tl , z im ) 
e,: random error of logpi ~ N( 0 , a 2 ). 

p estimated expected proportion defective of lot i at time t given 

(■2-il 1 Zim) 

t\\ estimated expected time that quality of lot breaches a predetermined 
level of proportion defective p given (.z,-!, z tm ) 

/ 3 0 , Pi', maximum likelihood estimators of /? 0 and / 3 t - 
logfi: an N xl matrix of logfii s, (logPi,..,logP^)' 

< 5 ,-: estimation error associated with logfii, < 5 , ~ N( 0 , if) 

7 q , rf: variance of fio and that of p t 

Tq, tJ: maximum likelihood estimators of and rf 

W: an NxN diagonal matrix with diagonal elements, l/(o r2 + t f/P?)'s 

7: maximum likelihood estimator for 7 

Mathematical formulations regarding the random effect logistic model 
are as follows: 

In the within- individual model, proportion defective of lot i at the jth 
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inspection, p^ is assumed to follow a logistic function of time t^: 

For i = 1, N, and j — 1 , n,- 

Ptj — exp(j3 o + (3itq)/(l + exp(j3o + Pitij)) (1) 

or 

log(Pij/(l - pij)) =P 0 + Pitij ( 2 ) 

In the between-individual model, deterioration rate of ammunition lot 
i, Pi, is represented as a function of covariates (e.g., depot condition, and 
vendor information) as well as a random error (e,): 



Pi ~ f (7l 2 il + •• + 7m z im + ^i) (3) 

where e, follows independent N(0,cr*). The proportion defective of ammu- 
nition lot would not decrease as time goes on. The following function is 
used to ensure positive deterioration rates, i.e. , positive Pi s: 

Pi = exp(yiZi] + .. + 7 m z tm + e.) (4) 

or 

l\ Z i\ T •• T 7m^im T Ci (5) 

After replacing Pi in ( 2 ) with that in (4), one could estimate regression 
parameters (Po, an d cr 2 ), by maximizing the marginal likelihood 

of the data based on the combined model. However, exact solutions for the 
normal equations based on the combined model are analytically and com- 
putationally infeasible (Stiratelli et al. (1984)). Indicating these problems, 
Stiratelli et al. (1984) discussed alternative ways of estimating parameters 
(EM algorithm and a two stage estimation). When P, in ( 2 ) is replaced with 
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that in (4), random error e,’s appear as coefficients of the regression model. 
In EM algorithm they are treated as missing observations and are esti- 
mated by cycling back and forth between the E-step and the M-step until 
convergence is reached. When a two-stage estimation method is used, the 
estimation procedure for the within-individual model is separated from that 
for the between-individual model. For instance, in the within-individual re- 
gression model, deterioration rates of a set of individual ammunition lots 
are estimated for all experimental lots. In the between-individual regression 
model, the estimated deterioration rates are used as values of a dependent 
variable and they are regressed on depot characteristics and vendor infor- 
mation associated with each individual lot. 

In general, a two-stage estimation method is computationally less inten- 
sive and analytically easier than the EM algorithm (Stiratelli et al. (1984)). 
In the following section, a two-stage estimation is described in order to for- 
mulate a prediction model for ammunition deterioration. 
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4. A TWO-STAGE ESTIMATION 

First, in order to estimate deterioration rates in the within-individual 



model (1), the following likelihood function of x^k conditional on /? 0 , /?!,..,/? at 
is formulated: 

N n, k) N n, / v 



i = n n ri ppv - = n n ( f .Wo - <«> 



»=iy=ifc=i 



i=iy=l VH 



where is as in equation (1). 

Maximum likelihood (ML) estimates, /?o, /?i,--Av are obtained by solv- 
ing a set of normal equations of log likelihood function of (6) with respect 
to j3 0 , /?i,..,/?jv. Estimated parameters (/?,’s) are in turn assumed to be ap- 
proximately normal with means equal to the true individual parameters 
(/?,•’ s) and variances (r 2 ’ s). Diagonal elements of the inverse of the nega- 
tive information matrix evaluated at (/?o, /?i,..,/W) provide the estimated 
variances of (/? 0 , which are (fj, ff, t£). 

Secondly, resultant /?, replaces unobservable /?,• in the between-individual 
model. This replacement, however, brings the estimation error 6, to the 
equation (5): 



l°g(/?i) = 7l-2jl + •• + 7 mZim + (7) 

where <5, is statistically independent of e,. It is assumed that <5, follows in- 
dependent Normal distribution with mean 0 and variance if//? 2 (var (<?,•) = 
var (logfii \ fii) ~ var ((/?,- — /?,)//?,) = rf //??). The var(<f) is approximated by 
ff//? 2 and /og/?* is marginally distributed as independent N(7 1 2 ll + , ..,+7 m 2 :,- m , 
<t 2 + ff//?f). Standard normal theory can be used to estimate the between- 
individual model parameters qj’s and <r 2 which maximize the following log 
likelihood function at the estimated rf : 
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N N ^ m 

LL = -J2 l og(v 2 + Ti)~ ~ Y1 'YiZa) 2 /(v 2 + *? /P 2 ). (8) 

i=l i=l 1= 1 

The resultant MLE of 7 is 

7 =(Z'WZ)- 1 (Z'W(logf3)) (9) 

and the estimated variance of 7 is (Z'WZ) -1 where an NxN diagonal 
matrix, W, consists of diagonal elements, l/(< 7 2 + t 2 //3 2 )’s. In order to 
be able to estimate 7 , the number of lots (N) employed in the experiment 
must be lager than the number of covariates (m) to be used in the between- 
individual model. In addition, as N increases, accuracy related to the 
inference concerning 7 would increase. 

When these 7 /s replace tj's in (7), logPi can be predicted in terms of 

log(4i) = 7l Z il + •• + ( 10 ) 

Use of (10) to quality control of ammunition, vendor control and depot 
maintenance plan is illustrated in the next section along with a numerical 
example. 
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5. IMPLIMENTATION 

Unfortunately, actual data related to type, quantity, and depot condi- 
tion of ammunition were omitted from the GAO report (The Comptroller 
General (1982)) to keep it unclassified. In the Norwegian study cited in the 
earlier part of this paper, repeated measurements on the deterioration of 
ammunition was not exposed for all the experimental lots. In fact, only the 
summary statistics for deterioration of several groups of ammunition were 
available, respectively. In order to illustrate implementation procedures for 
the methods suggested in this paper, a numerical example is given based on 
the simulated data generated partly based on the parameters used in the 
Norwegian study (Eriksen and Stromsoe (1974)) and the guideline of the 
U.S. Army Ammunition Surveillance Procedure (Supply Bulletin SB 742-1 
(1988)). 

Table 1 contains information regarding the series of the number of defec- 
tive items ( Yjj) found in 20 ammunition lots ( N ). These lots are classified 
into four groups depending upon the vendor who provided ammunition and 
depot in which ammunition is stored: group 1 consists of 5 lots (i — 1, .., 5) 
which are provided by vendor A and stored in depot 1; group 2 consists of 
5 lots (i = 5, .., 10) which are provided by vendor B and stored in depot 1; 
group 3 consists of 5 lots (i =11, .., 15) which are provided by vendor A 
and stored in depot 2; and group 4 consists of 5 lots (i — 16, ..,20) which 
are provided by vendor B and stored in depot 2. All the lots are inspected 
annually based on sample size of 20 ( Y t] ) . Total number of inspections of 
each lot (n,) varies over 4 different groups: n, = 7 for i — 1, .., 5; n, = 8 for 
i = 6, .., 10; n, = 9 for i — 11, .., 15; and n, = 9 for i = 16, .., 20. 

Our goals are first to find the relationship between deterioration rates 
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and concurrent characteristics of ammunition lot (depot condition and ven- 
dor information); secondly to determine if any group of ammunition causes 
significantly different deterioration rates from the others; and finally to pre- 
dict the time for renovation of ammunition in terms of characteristics of 
depot and vendor. 

In order to achieve the first goal, a two-stage estimation is employed. 
First of all, the ML estimates for parameters in the within-individual model 
0o, Ply , $20 and Tq, t 2 are obtained based on ( 6 ). For this step, 
PROC LOGISTIC of a statistical package SAS (SAS (1989)) is used by 
applying the values of Xg, 1^-, tij in Table 2 to model (2). Resultant /3's 
and rf's are summarized in Table 2 while (3 0 turns out to be -4.9034 which 
indicates that the estimated initial proportion defective of ammunition, 
(exp(J3 o)/(l + exp(j3 o))), is about 0.74% regardless of sources of vendor. 
Sample patterns of the actual deterioration (ActuahAf,/ Yij verses t^) are 
given in Figure 1 along with those of the estimated deterioration (Fitted: 
(exp0 o + /3,ty)/(l + exp0 o +Pitji)) verses t t] ). 

The between-individual model is formed to relate deterioration rates to 
the ammunition characteristics, logpi is used as dependent variable and the 
four dummy variables (z^,--, 2 *) representing 4 combinations of depot and 
vendor characteristics are used as covariates without an intercept: z,i = 1 
for ammunition lots which belong to group 1 , z,i = 0 otherwise; zp_ = 1 for 
ammunition lots in group 2 , z a = 0 otherwise; z& = 1 for ammunition lots 
in group 3, z^ — 0 otherwise; z , 4 = 1 for ammunition lots in group 4, 2,4 = 0 
otherwise. Given such ( 2 , 1 , .., Z&) as well as (/ 5 -, and f?'s), an optimization 
package GAMS (GAMS, 1988) is applied to ( 8 ) in order to obtain 71 , .., 74 
and a 2 . Estimates are summarized in Table 2 along with standard errors of 
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7 i , 74 obtained from the square root of the diagonal elements of ( Z ' XV Z ) 1 

/ 1 5 0 0 0 \ 



where a 20x4 covariate matrix Z = 



1 5 0 0 
0 1 5 0 



; 1 5 =(1,1, 1,1,1)’ is 



0 
0 

0 0 0 1 5 / 

a 5 x 1 unit vector; and W is a 20x20 diagonal matrix consisting of diagonal 
elements l/(<x 2 + r, 2 //? 2 ). Values of 77 and j3 }' s are given in Table 1 while 
<x 2 =0.012. 



Estimated 72 (/ = 1 , ..,4) is the weighted average of log(3 ' { s which belong 
to the group /. A large 72 indicates fast deterioration. For instance, the 
estimated average log deterioration rate of ammunition lots in group 2 (- 
0.8336) is lager than that of ammunition lots in group 1 (-0.9022). To test 
if any group characteristic causes significantly faster deterioration rate than 
another, a (1 — a)100% confidence interval for Z '7 can be used: 



Z' t(a/ 2, N - p)\[v(Z' WZ)^~1 



( 11 ) 



where Z is a 4x1 vector which contrasts a pair of two group characteristics; 
t is a student t distribution with degrees of freedom N — p. 

In order to test if the difference in two expected log deterioration rates 
between group 1 and group 2 is significant, one can use V = ( 1 , — 1 , 0 , 0 ) to 
obtain a 95% confidence interval for 71 — 72 : 



(1, -1, 0, 0)(— 0.9022, -0.8336, -0.7184, -0.6567/ 



±f(0.025, 20 - 4)^(1, -l,0,0)(Z'lVZ)- 1 (l, -1,0,0)' 

0.00708 0 0 0 



where ( Z ' WZ)~ X = 



\ 



0 0.00473 0 0 

0 0 0.00278 0 

0 0 0 0.00226/ 
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The resultant interval (—0.0686 ± 2.12(0.10870)) contains 0 which im- 
plies that at a = 5% vendor effect is not significant on the deterioration of 
ammunition stored in depot 1. Similar comparisons can be made by mod- 
ifying I's. Complete analyses of 95% confidence intervals for other pairs 
are summarized in Table 3. Overall, significant differences are not found 
except for the difference between group 1 and group 4. 

The resultant estimates from (10) can also be applied to the analyses of 
ammunition lots which are not used in the experiment. In order to predict 
the expected quality of a new ammunition lot P at time t in terms of depot 
and vendor characteristics associated with lot P the following p ) is used: 

Pi> — ex P0o + /^i’f)/(l + exp(f3 o + /?,-'<)) (12) 

where — exp(')iz i 'i +, , +7m2, 7 m)- For instance, the predicted proportion 
defective of ammunition lot P provided by vendor A after 5-year storage in 
depot 1 is 

p \ = exp(— 4.9034+ea:p(— 0.9022)x5)/(l+ea:p(— 4.9034+ea:p(— 0.9022)x5)) 

=0.0534. (13) 

In addition, one can estimate the expected time t) when the quality of 
ammunition lot { reaches a predetermined level p: 

t) = [log(p/( 1 - p)) - Po]/exp(y i Za +, , +n/ m z im ). (14) 

Expected time t\, when the quality of ammunition lot P provided by 
vendor A and stored in depot 1 reaches p=0.5 is estimated as 

t) = [log( 0.5/(l - 0.5)) + 4. 9034]/ea:p(— 0.9022) = 12.1. (15) 
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For the four groups considered in this example, p ) when t= 5 and t\ for 
p = 0.5 are given in Table 4. 
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6. SUMMARY AND DISCUSSION 

A random effect logistic regression model was used to formulate a predic- 
tion model for ammunition deterioration rate in terms of depot and vendor 
information. It was assumed that once lots are selected for initial inspec- 
tion, they will be inspected repeatedly over a certain time period without 
any rectification. The proportion defective was considered a measure of 
ammunition deterioration. Regardless of different sources of vendors, the 
initial quality of ammunition was assumed to be homogeneous as was its 
eventual deterioration with the passage of time. The log odds ratio due to 
unit time passage was defined as a deterioration rate which determines the 
shape of the logistic curve. Deterioration rates of different lots of ammuni- 
tion often vary and their variation was assumed to be explained partly by 
a function of depot characteristics and vendor information associated with 
each lot. 

A random effect logistic regression model was used to predict the dete- 
rioration rate of ammunition in terms of concurrent characteristics of the 
ammunition lot. A two-stage analysis was employed to estimate deterio- 
ration rates and relate them to the characteristics involved in ammunition 
lot. Using the approach of this paper, one can determine the appropri- 
ate time for reorder or renovation of ammunition before the quality reaches 
unacceptable levels when depot and vendor information associated with am- 
munition lot is available. By comparing the expected deterioration rates 
of ammunition stored in different locations one can improve depot mainte- 
nance policy. The similar approach can be used for vendor control. 

In the two-stage model used in this paper, we considered a certain type 
of ammunition and its deterioration was compared in terms of vendor and 
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depot characteristics. In order to compare deterioration rates of different 
types of ammunition, the following approach can be used: First, use a group 
of relatively homogeneous types of ammunition in the within-individual 
model, and secondly, classify them not only in terms of depot condition 
and vendor information but also by types of ammunition. 

In order to increase the accuracy of resultant estimators of a two-stage 
estimation, sufficiently large number of inspection which captures the no- 
ticeable deterioration of ammunition would be necessary over many exper- 
imental lots. However, such data may not be readily available in reality. 
When a sufficient number of repeated measurements of actual data is not 
available, accelerated testing may be an alternative means of meeting the 
time and budgetary restrictions of the experiment. For instance, envi- 
ronmental tests under simulated climatic extremes are an integral part of 
normal engineering tests. These tests include accelerated depot conditions 
at the test site when the characteristics of the test item necessitate the 
gathering of such data. Nonetheless, when only a small number of data 
is available, resampling procedure such as bootstrap method (Dalai et al 
( 1990 )) can be used to make inferences concerning deterioration rates at 
the cost of heavy computation. 

The measure of rates of ammunition deterioration considered in this 
paper was the proportion of defective items. When other variables such as 
muzzle velocity or water content are used as a measure of deterioration, 
modifications of the within-individual model may be necessary to take into 
account the different nature of the measure. However, the basic concept of 
a two-stage analysis still applies to any modified models. 

A restriction used in the between-individual model was that a lot of 
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ammunition would be stored in the same depot in the experimental period. 
A prediction model that reflects changes in depot conditions associated with 
a lot of ammunition is under investigation and is left for further research. 



20 



REFERENCES 

Brooke, A.; Kendrick, D. and Meeraus, A. (1988). GAMS: A User’s 
Guide , The Scientific Press, Redwood City, CA. 

The Comptroller General (1982). "DOD has Serious Problems with 
Care and Maintenance of Conventional Ammunition". AD A1 10764, Re- 
port to the Chairman, Subcommittee on Defense, Committee on Appropri- 
ations, House of Representatives. United States General Accounting Office, 
Washington D. C. 

Dalai, S. R., Fowlkes, E. B. and Hoadley, B. (1990). "Risk Analysis 
of the Space Shuttle: Pre-Challenger Prediction of Failure". Journal of 
American Statistical Association , 84, No. 408, pp. 945-967. 

Eckman, D. R. (1974). "Economic Comparison of Wood-Preservative 
Treated and Untreated 105mm Ammunition Boxes". AD A001532, Army 
Armament Command, Rock Island, IL. 

Eriksen S. and Stromsoe E. (1978). "An Examination of Deterioration 
of Ammunition by Storage, Norwegian Defense Research Establishment". 
AD A055897, Norwegian Defence Research Establishment, Norway. 

Korn, E. and Whittemore, A. (1979). "Methods for Analyzing Panel 
Studies of Acute Health Effects". Biometrics 35, pp. 795-802. 

Rao, C. R. (1965). "The Theory of Least Squares when the Parame- 
ters are Stochastic and its Application to the Analysis of Growth Curves". 
Biometrika 52, pp. 447-458. 

SAS (1990). SAS/STAT User’s Guide , Vol. 2, Version 6. Fourth Ed. 

Sohn, S. Y. and Mazumdar, M. (1991). "A Monte Carlo Study of Vari- 
able Selection and Inferences in a Two Stage Random Coefficient Linear 
Regression Model with Unbalanced Data", to appear in Communications 



21 



in Statistics, Theory and Methods A20. 

Stiratelli, R., Laird, N. and Ware, J. H., (1984) "Random-Effects Models 
for Serial Observations with Binary Response". Biometrics 40, pp. 961- 
971. 

Supply Bulletin SB 742-1, Inspection of Supplies and Equipment: US 
Army Ammunition Surveillance Procedures, Headquarters, Department of 
the Army, Washington DC. 

Vonesh, E. F. and Carter, R. L. (1987). "Efficient Inference for Random- 
coefficient Growth Curve Models with Unbalanced Data". Biometrics 43. 
pp. 617-628. 

Wilson, M. A. and Bird, R. (1986). "Defect Investigation of Cartridge 
Illuminating l-l/2inch". AD A173717, Technical Note, Department of De- 
fense, Defense Science and Technology Organization Materials Research 
Laboratories, Melbourne, Victoria, Australia. 



22 



ACKNOWLEDGMENT 

I am grateful to Mr. Larry Massa from Naval Weapon Support Center, 
Professors Judith Lind, Keebum Kang and Pat Jacobs from the Naval Post- 
graduate School for providing me with a great deal of practical information 
regarding ammunition deterioration and computational methodology. I al- 
so thank Ms. Hania Laborn for her editorial assistance. This research was 
partially sponsored by U.S. Naval Weapon Systems Center, Crane IN. The 
opinions expressed are solely those of the author. 



23 



Proportion Defective 




Actual 

o Lot 1 
□ Lot 6 
A Lot 11 
x Lot 16 



Fitted 




Lot 1 
Lot 6 
Lot 11 
Lot 16 



Year 



Figure 1. Deteriorating Patterns of Ammunition Quality 



Table 1: Number of Defectives X i} Observed in Sample Size of 20 at the 
jth Inspection of Lot i 
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Table 2: Fitting the Between-Individual Model 
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Table 3: 95% Confidence Intervals for the Difference between Expected Log 
Deterioration Rates of Two Groups 
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Table 4: Predicted Proportion Defective at year 5 and Estimated Time for 
p = 0.5 
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